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The leading semiclassical estimates of the electromagnetic Casimir stresses on a spherical and a 
cylindrical metallic shell are within 1% of the field theoretical values. The electromagnetic Casimir 
energy for both geometries is given by two decoupled massless scalars that satisfy conformally 
covariant boundary conditions. Surface contributions vanish for smooth metallic boundaries and 
the finite electromagnetic Casimir energy in leading semiclassical approximation is due to quadratic 
fluctuations about periodic rays in the interior of the cavity only. Semiclassically the non- vanishing 
' Casimir energy of a metallic cylindrical shell is almost entirely due to Fresnel diffraction. 

^ . 

■ 

^ ■ I. INTRODUCTION 

^ ' Casimir obtained his now famous attractive force between two neutral metallic plates [ij by considering the boundar" 
. conditions these impose on the electromagnetic field. Half a century later, his prediction was verified experimentally 
\^ ' to better than 1%. 

^-H , Twenty years after Casimir's work, Boyer calculated the zero-point energy of an ideal conducting spherical shellQ. 

Contrary to intuition derived from the attraction between two parallel plates, the sphere tends to expand. Boyer's 
result has since been improved in accuracy and verified using a number of field theoretic methods 0-01 ~ even though 
O ^ there may be little hope of observing this effect experimentally in the near future Q. 
^ I ^ ' Since field theoretic methods require explicit or implicit knowledge of cavity frequencies, they have predominantly 
, been used to obtain the Casimir energies of classically integrable systems. Thus, in addition to a spherical cavity, the 
d ' electromagnetic Casimir energies of dielectric slabs metallic parallelepipeds [T2|-[l^ and long cylinders 
^ , have been computed in this manner. 

1 ^ , '• However, most systems are not integrable and often cannot even be approximated by integrable systems. It thus 

', is desirable to develop reliable methods for estimating Casimir energies of classically non-integrable and even chaotic 
T-H ■ systems. Balian, Bloch and Duplantier calculate Casimir energies using a multiple scattering approximation to the 
^ ' Green's function(2P, 21]. This approach does not require knowledge of the quantum mechanical spectrum and the 
geometric expansion in principle is exact for sufficiently smooth and ideally metallic cavities. However, except for some 
' integrable systems, it in practice is often difficult to carry the multiple scattering expansion beyond its first few terms. 
. The relative importance of terms in this expansion also is difficult to assess a priori. In [2^ a semiclassical method was 
proposed to estimate (finite) Casimir energies. It is based on Gutzwiller's trace formula[23] for the response function 
and is suitable for Casimir energies of hyperbolic and chaotic s vstems [2^ [25| with isolated classical periodic orbits. 

Although not exact in general, the semiclassical approximation associates the finite (Casimir) part of the vacuum 
energy with optical properties of the system. It captures aspects of Casimir energies that have been puzzling for some 
time [26]. Path integral methods |27h31| in principle allow one to obtain Casimir forces to arbitrary precision. Due 
to unresolved renormalization problems, these methods have not yet been used to study the self-stress on cavities. 
We here use semiclassical methods adapted to classically integrable systems to estimate and analyze the Casimir 
^ ' self-stress of a spherical and a cylindrical shell. 
5^ , The simplicity, transparency and surprising accuracy of this approximation is first demonstrated on Boyer's 
problem [J-Q, the electromagnetic Casimir energy of a spherical cavity with an (ideal) metallic boundary. The semi- 
classical analysis of this problem is an order of magnitude simpler than any given previously and the positive sign 
is related to caustic surfaces of second order. In sect. 4 the semiclassical estimate for the electromagnetic Casimir 
self-energy of a perfectly conducting cylindrical shell is reexamined by converting the sum of WKB-estimates for the 
eigenfrequencies to the dual sum over periodic orbits. It was previously [H, HI] found that the Casimir self-energy of 
a cylindrical metallic shell vanishes to leading semiclassical order. It in fact vanishes only because the upper bound 
of a particular Fresnel integral is ignored in this approximation. Demanding this physical bound in the longitudi- 
nal momentum fraction, the semiclassical estimate also is rather accurate for the self-stress of a metallic cylindrical 
shell and within |% of its field-theoretic valuefTsI]. The sign of the self-stress is again determined by the presence 
of caustics and the associated Maslov-Keller indices 34, 35]. However, contrary to the spherical case, this self-stress 
semiclassically is primarily due to Fresnel diffraction effects. A discussion of the results and a critical assessment of 
difficulties remaining for a semiclassical interpretation of some Casimir self-energies concludes this article. 
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II. THE DUAL PICTURE: ELECTROMAGNETIC CASIMIR ENERGIES OF INTEGRABLE SYSTEMS 

AND PERIODIC RAYS 

Integrable systems may be semiclassically quantized in terms of periodic paths on invariant torifslj - in much the 
same manner as Bohr first quantized the hydrogen atom. Ahhough in general not an exact transformation, classical 
periodic orbits on the invariant tori are dual to the mode frequencies in a semiclassical sense. Applying Poisson's 
summation formula, the semiclassical Casimir energy (SCE) due to a massless scalar may be written in terms of 
classical periodic orbits [2 3t i26> ,37, . .38.] , 



The components of the d-dimensional vector I in Eq.([T]) are the actions of a set of properly normalized action- angle 
variables that describe the integrable system. The exponent of the integrand in Eq.([T]) is the classical action (in 
units of h) of a periodic orbit that winds nii times about the i-th cycle of the invariant torus. H{T) is the associated 
classical energy and /3rn is the Keller-Maslov index [s^ [ssj of a class of periodic orbits identified by m. The latter is a 
topological quantity that does not depend on the actions I. To leading semiclassical order, the (primed) sum extends 
only over those sectors m with non-trivial stationary points (the classical periodic paths of finite action) (see below) . 
The correspondence in Eq.([T]) can only be argued semiclassically [H, 113, HI] and the integrals on the RHS should 
therefore be evaluated in stationary phase approximation (sp) only. 

The semiclassical spectrum of a massless scalar is exact for a number of manifolds without boundary [s^ and the 
definition of the semiclassical Casimir energy (SCE) by the RHS of Eq.(IT]) coincides with the Casimir energy of 
zeta-function regularization in these cases. It also is exact for massless scalar fields satisfying periodic-, Neumann- or 
Dirichlet- boundary conditions on parallelepipeds [l^, [13, as well as for some tessellations of spheres [2^ l4l| . 
To physically interpret the finite SCE of a system, one has to consider the implicit subtractions in the spectral 
densityHH, 0,143]. 

Semiclassical contributions to the spectral density arise due to small fiuctuations about closed classical paths. These 
either are contractible to a point or not. in this dual picture local UV divergences are associated with fluctuations about 
arbitrary short contractible classical paths. If all local UV divergences can be subtracted unambiguo usly[2i|4l,|4l],the 
dependence of the remaining finite Casimir energy on macroscopic deformations of the system arises from fluctuations 
about classical closed paths of finite classical action only. 

In disjoint systems, the Casimir interaction at small separations usually is dominated by fluctuations about periodic 
orbits[26]. If there are no stationary classic al p eriodic orbits, the leading interaction typically is due to periodic orbits 
of extremal- rather than stationary- length[4J]. These correspond to diffractive effects that lead to relatively weak but 
sometimes rather interesting interactions. Examples of systems without stationary periodic orbits include the Casimir 
pendulum of[45j, perpendicular plates [29|. a wedge above a plate and, most recently, a rotational ellipsoid above a 
plate with a hole[46|. Diffractive contributions also become important when the separation between two systems is 
comparable to or larger than the smaller system. The original Casimir- Polder interaction between two atoms [4?!] or 
between an atom and a metallic plate ^48] fall in this category. 



But periodic (stationary or extremal) orbits need not dominate the Casimir energy when some non-periodic closed 
orbits of finite action are much shorter than any periodic ones. The Casimir force due to a massless scalar field 
satisfying Dirichlet boundary conditions on the plate of a hemispherical Casimir piston j49| is an example. Periodic 
orbits similarly may not give the largest contribution to the self-stress of a cavity for scalar fields satisfying Dirichlet 
or Neumann conditions. For sufficiently smooth boundaries like the spherical and cylindrical (toroidal) shell we are 
interested in, only arbitrary short classical paths that reflect just once off the boundary potentially can lead to UV- 
divergences. The divergence generally will depend on the local curvature of the boundary [43]. Even if this divergence 
is not logarithmic and can be subtracted unambiguously, the finite remainder of such short closed orbits can contribute 
significantly to the Casimir energy [4^. 

The following dimensional argument suggests that in a dimensionless regularization scheme the local ultraviolet 
divergence due to scalar fields satisfying Dirichlet or Neumann conditions on both sides of a smooth and infinitesimally 
thin even-dimensional surface in fact vanishes: barring other scales, the local contribution from a small {d — 1)- 
dimensional surface element dA to the divergence for dimensional reasons is of the form hcfe{Ri/ Rj) dA/R'^, where 
R is the principal (local) radius of curvature of the surface and is a dimensionless function of the (dimensionless) 
regularization parameter e and of ratios of the local curvatures only. The interior and exterior radii of curvature at 
the same point on the infinitesimally thin surface are of equal magnitude but of opposite sign. Local divergent surface 




n 



UV subtractions 



n 




(1) 



103 



contributions from the interior and exterior of the infinitesimally thin surface [with the same scaleless boundary 
conditions on both of its sides] thus cancel precisely for odd dimensions d. For a spherical surface this cancelation 
has been explicitly observed injs^l- In d = 3, finite Casimir energies have also been calculated for scalar fields and an 
infinitesimally thin cylindrical shell fl6l [Sll. [52j. The argument above suggests that surface divergences in fact cancel 
locally for any infinitesimally thin (and sufficiently smooth) even-dimensional boundary, regardless of its shape or 
whether Dirichlet or Neumann conditions are imposed. 

However, for scalar fields this cancelation occurs only for vanishingly thin even-dimensional smooth boundaries and 
in a dimensionless regularization scheme. In a physical cutoff scheme this only implies the absence of logarithmic 
divergences, and surface divergences do occur even for ideal boundaries. They can be unambiguously subtracted fssl- 
[56j or, equivalently, absorbed in parameters describing physical properties of the surface. Although the subtractions 
are not ambiguous for even-dimensional ideal interfaces, the remaining finite contribution to the Casimir energy due 
to closed non-periodic classical paths could be sig nificant and Eq.(IT|) in this case can give an inadequate estimate of 
the Casimir self-stress due to scalar fields [see [la El] and sect. -4 for some formally sub-leading semi-classical surface 
contributions that have been omitted in Eq.((T]).] 

Fortunately the finite electromagnetic self-stress of a smooth and infinitesimally thin perfectly conducting boundary 
is predominantly due to the stationary points of the integrand in Eq.([T]). The electromagnetic Casimir self-energy 
of a closed, smooth and perfectly metallic shell may be decomposed into the contributions from two massless scalar 
fields - one satisfying Dirichlet's, the other Neumann's boundary condition[69] on the surface [2l|. Semiclassically, 
these boundary conditions are enforced by a phase lag of 7r(0) for the Dirichlet (Neumann) scalar at each specular 
refiection of the classical path. Since both scalar fields satisfy the same wave equation, only fluctuations about 
classical trajectories with an even number of reflections contribute to the electromagnetic spectral density of an ideal 



metallic cavity [21[. Potentially divergent contributions of the Neumann- and Dirichlet- scalar that arise from closed 
contractible classical paths that reflect just once off the metallic boundary in this case cancel each other exactly. The 
length of any closed classical ray with an even number of reflections off a sufficiently smooth cavity (without sharp 
corners) thus is bounded below by geometrical characteristics of the cavity - such as its minimal radius of curvature. 
The closed classical rays of minimal length with an even number of reflections off such cavities are periodic and 
correspond to either a stationary point or an extremum of the classical action. For the cylindrical- and spherical- 
shells we will be considering, extremal periodic ray s creep about the exterior of the cavity a number of times. They 
lead to rather small diffractive corrections [s^ 13 that will be ignored. The stationary points of the action in Eq.([J) 
give the main contribution in the electromagnetic case. They correspond to periodic trajectories inside the cavity 
characterized by their winding number and (even) number of reflections such as those shown in Fig. la. 

The previous argument applies equally well to the contribution of any pair of decoupled scalar fields satisfying 
Dirichlet and Neumann boundary conditions and it may not be apparent why the SCE should give a particularly good 
approximation for the self-stress of a metallic shell due to the electromagnetic field. The exact field theoretic Casimir 
stress caused by two massless scalar fields solving the same Helmholtz equation as the transverse electromagnetic fields 
but satisfying Neumann- (instead of Robin-) and Dirichlet- boundary conditions on the spherical shell not only differs 
in sign but also is an order of magnitude larger than the electromagnetic one'TG','!^,!!^: —0.220967 . . . hc/R instead of 
the -f 0.046 17 ... ftc/i? in the electromagnetic case^St ifij]. The reason for this difference is that the Neumann boundary 
condition on a spherical shell is not conformally covariant [STj for a scalar satisfying the Helmholtz equation whereas 
the Robin boundary condition is[70j. On an intrinsically flat boundary like a cylindrical shell, the electromagnetic 
Casimir energy of a metallic shell indeed decomposes into contributions from two scalar fields satisfying Dirichlet and 
Neumann boundary conditions[l^. The classical action of a massless scalar particle is conformally invariant and the 
semi-classical approximation indeed reproduces [2^ the Casimir energy of a conformal scalar field on curved manifolds 
without boundary such as that of a three-dimensional sphere S3. Since specular reflection and phase lag do not 
depend on the curvature of a surface, we conjecture that for a smooth but curved boundary the SCE approximates 
the Casimir energy of scalars satisfying conformally covariant Dirichlet and Neumann conditions on the boundary. 
For boundaries with non- vanishing curvature, the latter correspond to Robin-like conditions. 



III. SELF-STRESS OF A SPHERICAL METALLIC SHELL 



Classically a massless particle in a spherical cavity is an integrable system, but the semiclassical spectrum is only 
asymptotically correct. One therefore cannot expect the semiclassical approximation to be exact in this case. It 
nevertheless will be surprisingly accurate. The SCE is obtained by performing the integrals of Eq.([T]) in stationary 
phase and has a very transparent interpretation in terms of periodic orbits within the cavity only. The sign of the SCE 
of a spherical cavity in particular will be quite trivially established and the discrepancy of 1% with the fleld-theoretic 
results may very well largely be due to diffractive corrections from creeping orbits that wind about the exterior of the 
sphere. As argued above there are no (potentially divergent) local contribution to the Casimir energy from such an 
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idealized surface in the electromagnetic case - its local surface tension in fact vanishes [21], '60'|. The only subtraction 
in the spectral density required for a finite Casimir energy with ideal metallic boundary conditions is the Weyl 
contribution proportional to the volume of the sphere. This subtraction corresponds to ignoring the m = (0, 0, 0) 
term in the sum of Eq.(IT]). The remaining difficulty in calculating the SCE is a convenient choice of action-angle 
variables. For a massless scalar in three dimensions satisfying boundary conditions with spherical symmetry, an 
obvious set of actions is the magnitude of angular momentum, I2 = L, one of the components of angular momentum 
^3 = Lz and an action Ii associated with the radial degree of freedom. 

Since the azimuthal angle of any classical orbit is constant, the energy E ~ H{Ii,l2) of a massless particle in a 
spherical cavity of radius R does not depend on I3, = L^- In terms of this choice of actions, the classical energy is 
implicitly given by, 

./.+/.arccos(^).^yC(^. (2) 

The branches of the square root and inverse cosine in Eq.([2]) are chosen so that Ii is positive. It is convenient to 
introduce dimensionless variables 



A = 2ER/{hc) and z = ch/iER), (3) 

for the total energy (in units of hc/{2R)) and the angular momentum (in units of ER/c) of an orbit. Note that 
z G [0, 1] and that the semiclassical regime formally corresponds to A ^ 1, i.e. to wavelengths that are much less than 
the dimensions of the cavity. Using Eq.Q and the definitions in Eq.Q, the angular frequency of the radial motion 
is = (dE/dh)-^ = i?Vl - ichl{ER)Y/{^c) = {R/ttc)VT^. 

With the help of Eq.® and the definitions of Eq.([2]), the semiclassical expression in Eq.([l]) for the Casimir energy 
of a massless scalar field satisfying Neumann or Dirichlet boundary conditions on a spherical surface becomes. 
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The integral over I3 has here been performed in stationary phase approximation. Because the Hamiltonian does not 
depend on I^, only periodic orbits with — contribute [26] significantly in stationary phase. Since —I2 < ^3 < ^2, 
one has that J dl^ ~ 2/2 = Xz. The factor 2/2 accounts for the 2{l + l/2)-degeneracy of a state with angular 
momentum L = I + 1/2 = 12- By taking (4 times) the real part in Eq.(j4]) one can restrict the summations to non- 
negative integers and choose the positive branch of the square root- and inverse cosine- functions in the exponent. 
The primed sum here means that the summand is weighted by half if one of the integers vanishes and the absence of 
the n = w = term. The Keller-Maslov index (3{n, w) of a classical sector depends on whether Neumann or Dirichlet 
boundary conditions are satisfied on the spherical shell and is obtained below. 

For positive integers w and n, the phase of the integrand in Eq.Q is stationary at z = z{n, w) G [0, 1] where, 

= — n arccos(z) -t- WTT 

z{n, w) — cos{wTr/n), n > 2'w > 1 . (5) 

Restrictions on the values of w and n arise because arccos(2) G [0, tt/2] on the chosen branch. The phase is stationary 
at classically allowed points only for sectors with n > 2^ > 1. Semiclassical contributions to the integrals of other 
sectors arise due to the endpoints of the z-integration at 2: = and z — 1 only. These diffractive contributions are 
of sub-leading order in an asymptotic expansion of the spectral density for large A and will be ignored. Note that 
w ^ w + n just amounts to choosing another branch of the inverse cosine. 
The classical action in sectors with stationary points is, 

Sciin^w) = hXn sin{wTr / n) (6) 
= {E/c)2nRsiii{wTT/n) = (E/c)L{n,w) , 

where L(n, w) is the total length of the classical orbit. Some of these classical periodic orbits are shown in Fig. 1. The 
integer w in Eq.® gives the number of times an orbit circles or winds about the origin. The integer n > 1 gives the 
number of reflections off the spherical surface (windings of the radial motion). As indicated in Fig. 1, the envelope of 



105 



the set of classical periodic orbits in the (n, t(;)-sector is a caustic surface and a double covering is required for a unique 
phase-space description ^3]. The two sheets are joined at the inner caustic [indicated by a dashed circle in Fig. 1] and 
at the outer spherical shell of radius R. Every orbit that passes the spherical shell n times also passes the caustic n 
times. The cross-section of a bundle of rays is reduced to a point at the spherical caustic surface. The caustic thus 
is of second order and associated with a phase loss of tt every time it is crossed. At each specular reflection off the 
outer shell Dirichlet boundary conditions require an additional phase loss of tt whereas there is no phase change for 
Neumann boundary conditions. Altogether the Keller-Maslov index of the sector (n, w) depends on n only and is 
given by, 

a, \ / 0: ^OT Dirichlet b.c. 

^ ' ' 1 2n, for Neumann b.c. ' ^ ' 

As noted in sect. 2, the electromagnetic field satisfying (ideal) metallic boundary conditions on a spherical shell, may be 
viewed as two massless scalar fields of vanishing conformal dimension, one satisfying Dirichlet and the other Neumann 
boundary conditions [2l|. Due to the Keller-Maslov phases of Eq.([7]) only sectors {n,w) with even n = 2k > 2w > 1 
contribute j21,J to the SCE in the electromagnetic case. 




FIG. 1: (Color online) Classical periodic rays of a spherical and cylindrical cavity, a) The shortest primitive rays corre- 
sponding to sectors (n, w) € {(2, 1), (4, 1), (6, 2)} that contribute to the electromagnetic SCE. b) closed paths in sectors 
(n, w) G {(1, 0), (3, 1), (5, 1), (5, 2)} that reflect an odd number of times olT the surface and whose contribution to the elec- 
tromagnetic SCE vanishes. Caustic surfaces are indicated as thin circles. The dashed part of any trajectory is on one sheet and 
its solid part on the other of a two-sheeted covering space. The "phase space" of the (5, 2) sector is indicated by the hatched 
area. Note that the caustics are of 2"^* order for a spherical cavity but of order for a cylindrical one. 



The classical action vanishes for sectors with w — or n = and these sectors do not contribute to the SCE in 
leading approximation. The n = sector corresponds to paths that do not reflect off the shell and thus gives rise 
to volume contributions that are subtracted. Eq.© implies that extremal paths in the {n > 0,w — 0) sectors have 
maximal angular momentum 1 — z = lc/{ER). These are great circles that lie wholly within the spherical shell. 
These paths are not stationary and lead to diffractive contributions that we will ignore here. For (n > 0, w > 0) the 
curvature of the exponent at z(n, w) is finite. 



z arccos(2)) + wkz] 
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whereas it diverges in sectors with w ~ 0. The behavior of the exponent for z 1 in this case is. 



Vi-z^ — z arccos(z) 



(9) 



Quadratic fluctuations about the classical orbit with w — ^ thus are of vanishing width and these sectors do not 
contribute in stationary phase approximation. To leading semiclassical accuracy, the Casimir energy of a spherical 
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cavity with an ideal metallic boundary therefore is, 
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1 - 15V2 cos(f^) 
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- 0.04668...^. (10) 
-ft 

This semiclassical estimate is only about 1% larger than the best numerical valuej^ 0. 04617. ..ftc/i? for the electro- 
magnetic Casimir energy of a spherical cavity with an infinitcsimally thin metallic surface. Note that the contribution 
from the {211!, w) sectors had to be considered separately in Ea.([TU| since the measure dzz vanishes at the stationary 
point z{2'w,'w) — cos(7r/2) = of the integrand, which is an endpoint of the integration domain. As can be seen in 
Fig. la), the classical rays of (2i(;, z/;)-sectors go back and forth between antipodes of the cavity and pass through its 
center - they have angular momentum L = 0. 

The shortest primitive orbits give somewhat less than half [1 / (167r) ^ 0.02] of the total SCE of the spherical cavity - 
much less than the 92% they contribute to the Casimir energy of parallel plates. The reason is that contributions only 
drop off as 1/n^ rather than like l/n* as for parallel plates. The length of an orbit in the (4, l)-sector (the inscribed 
square in Fig. la) furthermore is just a factor of \/2 longer than an (2, l)-orbit [which in turn is a factor of l/\/2 
shorter than an (4, 2)-orbit]. To estimate the magnitude of the contribution of any particular sector one has to take 
the available phase space as well as the ray's length into account. Thus, although the length of a (2n, l)-orbit tends 
to 2ttR for n — )■ oo, the associated phase-space (essentially given by the volume of the shell between the boundary of 
the cavity and the inner caustic) decreases like 1/rt^. This accounts for the relatively slow convergence of the sum 
in Eq. fTUl) . To achieve an accuracy of 10^^, the first 50 terms of the sum were evaluated explicitly and the remaining 
contribution was estimated using Richardson's extrapolation method. However, note that the semiclassical 2-reflection 
coefficient is just 2/5-ths of the leading field-theoretic coefficient for a dilute dielectric-diamagnetic spherical shell j31^. 
We will encounter a similar discrepancy in the case of a cylindrical shell and will discuss it further in the conclusion. 



IV. SELF-STRESS OF A CYLINDRICAL METALLIC SHELL 



The example of a spherical cavity shows that one may obtain electromagnetic Casimir self-energies rather accurately 
by considering only small fluctuations about their classical periodic orbits. Since the classical periodic rays of finite 
length for a long cylindrical shell are the same as for a sphere, one would expect that a semiclassical calculation of 
the self-stress is just as straightforward for a cylindrical cavity. 

This is not the case. We verify below that the electromagnetic SCE of an ideal metallic cylindrical shell vanishes 
in stationary phase approximation fs^ [33| . The contribution from any periodic orbit to the electromagnetic SCE 
vanishes in this approximation for the same reason that it is positive for a spherical cavity - due to optical phases. 
However, going beyond stationary phase approximation and including certain Fresnel diffraction effects, we again 
obtain a very good approximation to the Casimir self-stress of a metallic cylinder from quadratic fluctuations about 
classical periodic paths. 

The cylinder appears particularly suited for a semiclassical analysis in terms of massless scalars because the trans- 
verse electric- and magnetic- modes satisfy Dirichlet and Neumann boundary conditions on the cylindrical surface. 
The classical system is integrable and we could directly employ the formalism of Berry and Tabor [26l [37j or an ex- 
tended Gutzwiller approachjH, to obtain the SCE in much the same way as we did for the sphere. However, 
to better understand why the SCE vanishes in this approximation and improve upon it, we first obtain the dual 
expression for the Casimir self-energy of a cylindrical cavity directly from the semiclassical (WKB) estimate of the 
eigenvalues of the scalar fields. The Casimir energy of an infinitcsimally thin metallic cylindrical shell of radius R 
within a much larger cylinder of fixed radius i?> ~ oo in principle is given by the i?-dependent part of the zero-point 
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energy, 



+ + i?>) - + 4(0, i?>)| , (11) 

J D.N 



where {Kn{R<, R>)}d,n is the spectrum of wave numbers for a scalar field satisfying Dirichlet(Z?) or Neumann(A^) 
boundary conditions on a 2-dimensional annulus with inner radius i?< and outer radius i?> . 

The asymptotic heat-kernel expansion[6^ implies that the subtracted expression in Eq. ljlll) is finite: the potentially 
logarithmic divergence proportional to the average of the 3rd power of the (extrinsic) curvature of the cylindrical 
surface is canceled. Milton and DeRaad calculated this finite self-energy exactly some time ago[l^. Since they 
uniformly approach the exact wave numbers sufficiently rapidly, the SCE obtained by replacing the exact eigenvalues 
in Eq. pH) by their WKB-estimates also is finite. For i?> ^ oo all periodic orbits in the annulus have a length of 
C'(i?>) and the only finite contribution to the SCE from the annulus is due to (diffractive) creeping orbits that wrap 
about the inner cylinder. We will neglect this small contribution to the SCE and consider only radius-dependent 
semiclassical contributions from the inner cylinder. The transverse wave numbers Kne{0,R) = Xnt/R of interior 
modes of the cylinder in WKB-approximation are positive solutions of [ssl. [63| . 

ft{x,,p)^^{n + \±\), forn = 0,l,...;^-0,l,... , (12) 

where the (-1-)- and (— )-sign corresponds to Dirichlet and Neumann boundary conditions respectively and, 

fi{x) = V2:2-^2_^arccos(^/a;) . (13) 

This semiclassical (Debye) approximation generally gives the zeros of Bessel functions of the first kind and their 
derivatives to better than 1%. There is but one notable exception: the zero of Jq{x) = Ji{x) at x = corresponds to 
a WKB value of xqo = 7r/4. Note that fi{x = €) = Q and all semiclassical wave numbers satisfy Xni > i- 

Inserting the approximation in Eg. dT^ and Eq. p^ for the eigenvalues, the i?-dependent contribution of interior 
modes to the SCE of a conducting cylindrical shell is the finite i?-dependent part of, 

n^ — oo i—Q 
oo noo 



X dy / dxy/y^+x^f'/.{x)e'^"'f'^'''> 

Jo J I 
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dy dxy^y^+x^dx\n[cos{2ff{x))] . 

(14) 



The summation over n has been performed in the last expression. It shows the equivalence of the present approach with 
one based on the generalized argument principle [gI] with the function, cos(2/£(x)) = ^{e~'^^''^'^^ + ze'-^*^'^')(e~*'^*('^) — 
ie^ft(^)'j^ whose zeros are the WKB estimates for the mode frequencies for each partial wave, and a contour that 
runs from x — od to x ~ £ just below the real axis and returns to a; = 00 just above it. Although divergent 
because it includes only the contribution from the interior, it is reassuring that the expressions in Eg. dTi]) are not 
logarithmically divergent. The divergences therefore can be subtracted unambiguously (and in fact are canceled by 
exterior contributions we are not considering) and a semiclassical evaluation is possible. 
One arrives at Eg. (1141) by applying Poisson's resummation formula, 

00 00 

n— — 00 n— — 00 

to the dimensionless (and scaled) semiclassical spectral densities, 

00 00 

Pd/n{x) = EE'^(^~^"'^) 

e=0 n=Q 
00 00 

^J2e{x~£)n{x) ^ S{Mx)-{n+^±\)7T) . 
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The Heaviside function 6{x — £) in Eq.([T6]) ensures that fi is real. Using Poisson's relation once more in the form, 



oo oo 
D — n 111 — ''0 



de e^'^™'5(^) , (17) 



£—0 u) — — oo ^ 

the sum over partial waves of cavity modes is converted to one over the winding number of classical paths. In terms 
of the dimensionless wave number A = ya^^+y^ = ER/ {he) and the longitudinal- and angular- momentum fractions 
a = y/\ and z = Ijx one obtains, 



. ^ oo oo ^oo /■! 



(18) 
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^2 — ^2g2iA\/l — Q^(4n(\/l-2'''-zarccos z)+WTrz) 




We here used that only terms with even n contribute in Eq. (jl4p . Only the principal branch of the inverse cosine is to 
be considered here and the integrals are formal in the sense that they are to be evaluated to leading (non-vanishing) 
order in stationary phase approximation only. The integrals are finite in this restricted sense and their divergent part 
has been implicitly subtracted. They otherwise diverge (in all sectors) due to the behavior of the integrand at large 
A near a ^ 1, i.e. for large longitudinal momentum fractions. As explained in sect. 2 we are assured that all surface 
divergences cancel for metallic boundary conditions on a cylinder. We are interested only in the finite contributions 
arising from quadratic fluctuations about periodic classical trajectories (stationary points of the integrand in the 
(n 7^ 0, ui 7^ 0) sectors). 

The ^5(0) term in Eq.p7|) gives a surface correction [ssj to Ea. ([T5)) . As mentioned in sec. 2, such surface 
contributions cancel in the electromagnetic case. One may explicitly verify this by subtracting the n = sector and 
noting that 
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The stationary points of the z-integral in Eq. (jl8p for w < n > at Znw = cos {■Kw/n) are the same as for a spherical 
cavity and correspond to planar periodic orbits like those of Fig. la. To quadratic order the fluctuations about the 
stationary point at a = give rise to Fresnel-like integrals of the type (for generic 7 > 0), 

1 

^^gi7r7'x/T^ _ g»^7M da e-"'^'"'/2 ^ (20) 
"'0 

= 6--^' (C(7)-z5(7))/7. 

Extending the upper bound of the fluctuation integral in Eq. p0|) to 00 and thereby replacing the cosine- and sine- 
Fresnel integrals of Eq. (f20|) by their mean value of 1/2 gives a vanishing value for the SCE of a metallic cylindrical 
shell js^ [33}. The leading non- vanishing contribution in this case arises from the finite upper bound of the fluctuation 
integral in Ea. (l20l) . The SCE of a metallic cylindrical shell in this semiclassical sense is almost entirely due to Fresnel 
diffraction effects. 

The semiclassical evaluation of the integrals in Eq.(fT8|) again is facilitated by noting that one can take (4 times) the 
real part of the (n > 0, w > 0)-sector contributions and that the spectral density is analytic in the first quadrant. The 
X — iS, integral in Eg. (1181) may thus equally well be performed along the positive imaginary axis. The contour is closed 
in the complex A-plane by a large quarter circle on which the integral, or better, a regularized version of it, vanishes 
sufhciently rapidly. The integrand is real on the imaginary energy axis. For Q < w < n the stationary points of the 
integrand in Eq. (jl8p are at z — cos ^ and a — 0. They correspond to planar periodic rays of a cylindrical cavity 
like those in Fig. la with winding number w and 2n vertices. Expanding to quadratic order about these stationary 
points, we first integrate out fluctuations in z. Rescaling ^ ^/(l - the remaining ^— and a— integrations may 

then be carried out independently. Note that the analytic continuation of the integrals is possible only if we assume 
an upper bound at a = 1 — e. Starting with the expression of Eq. (|2ip this procedure gives the following semiclassical 
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estimate of £^J^ (with c. 
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7tt{7tt^ ~ 240) fcL 



276480 
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= -0.013594.. 



hcL 



The SCE of a conducting cyhndrical shell thus differs by less than 0.25% from the exact field theoretic value |15| 
of —0.0135613 ... ?icL/i?^ if a certain kind of Fresnel diffraction is included. It was necessary to go beyond the 
formally leading semiclassical approximation because the latter violates causality and vanishes. Enforcing the bound 
on longitudinal momenta, the semiclassical approximation again closely reproduces the Casimir self-stress on a metallic 
cylindrical shell. 




X 

FIG. 2: Integrating over the longitudinal momentum fraction. The upper solid curve shows the exact integral in Eq. (|22p 
over the longitudinal momentum fraction a, the lower solid line is its semi-classical approximation to quadratic order in the 
fluctuations. The dashed curve gives the simple exponential approximation e~'^^^'''. 

To better understand the approximation made, compare the following approximations to the Fresnel-like a-integral 
in Eq.CT: 

dae-"^i^ = I (H^iitx) - h{x)) - e^^^/^ 

dae-^(-4)_,,^^y|erf(.y|) (22) 

The first line expresses the a-integral occurring in Eq. ljlSp for values of the energy on the positive imaginary axis 
in terms of Struve- and Bessel- functions. The function e^'^^/'* gives the best exponential fit for small x ^ 0. The 
second line of Ea.(|22p is the semiclassical approximation to second order in the fluctuations. It also gives a uniform 
approximation to the integral that decays exponentially for x ^ oo (but as e*- — x/2)/x). As shown in Fig. 2, both 
approximations reproduce the integral well for small values of a: - where the result is sizable - but cut down its slow 
decay (~ l/s;^) at large x ~ cxi. As noted before, the power-like decay of the original integral gives rise to divergent 
contributions to the Casimir energy in every {n,w) sector [tJI. Assuming (and heat-kernel considerations show this is 
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possible for an ideal metallic shell [53|, |65|) all divergences are subtracted unambiguously, the leading non- vanishing 
semi-classical approximation apparently is quite a reasonable estimate of the finite part. One in principle could 
improve on the representation of the integral by including higher orders in the expansion of \/\ — o? about a = 
- but these contributions would be inconsistent with the (leading order) WKB approximation to the eigenvalues. 
We can to some extent assess the sensitivity of the SCE to the precise manner in which the power-like tail is cut 
off by comparing with the exponential fit — 7rx/4) for the a-integral: in this approximation the Casimir energy 
for a metaUic cylinder becomes (Ttt^ - 2m)hcL j {2%%t^'^ ^B?) = -0.013533 . . . hcL/R^ - well within the error of the 
semiclassical estimate. 



V. DISCUSSION AND CONCLUSIONS 



We here obtained the electromagnetic Casimir self-stress of perfectly conducting spherical and cylindrical shells 
in semiclassical approximation. This approach reproduces the field theoretic values 0, Q to better than 1%. The 
semiclassical description by two massless scalars in general gives the electromagnetic Casimir-stress of a metallic 
spherical shell rather than the sum of the Casimir stresses due to scalar fields satisfying Neumann and Dirichlet 
boundary conditions. For a spherical shell the latter has the opposite sign and is an order of magnitude larger 

[MS 

[59| . The semiclassical description is inherently conformal and thus may be best suited to describe scalar fields that 
satisfy conformally covariant boundary conditions |57j. It was previously observed 26] that the dispersion for massless 
scalar particles leads to a semiclassical description of conformally coupled scalar fields on curved spaces such as 6*3. 
They also appear to couple conformally to curved boundaries. "Neumann" boundary conditions are semiclassically 
imposed by requiring no phase change and specular reflection. They do not change under conformal rescaling of the 
boundary and in this sense are conformally covariant. For spherical boundaries they describe (massless) scalar fields 
satisfying conformally covariant [57] Robin conditions. This conjecture is supported by the fact that the low- lying 
semiclassical eigenfrequencies implied [tI] by Eq.([TU]) are closer to those of two scalar fields satisfying Dirichlet- and 
Robin- than Dirichlet- and Neumann- conditions. Robin boundary conditions approach Neumann boundary conditions 
at high frequencies and both would be similarly implemented in geometrical optics (i.e. with specular reflection and 
no phase loss). 

On a cylindrical surface in flat space ordinary Neumann and Dirichlet conditions are already conformally covariant. 
We estimated the electromagnetic Casimir self-stress of a metallic cylindrical shell by explicitly summing WKB- 
approximations to the eigenfrequencies of the two scalars in the dual picture using Poisson's resummation formulae. 
The dual of the principal quantum number gives the number of interactions with the shell. Summing over these one 
arrives at the general argument principle (see Ea. (jlip ) often used[63| to evaluate Casimir energies when the spectrum 
is only implicitly known. The difference here is that the roots of the characteristic function are WKB estimates of the 
eigenfrequencies. The dual to the sum over partial waves is the sum over windings about the origin. The stationary 
points of minimal action in each sector are classical periodic trajectories. We evaluate the integrals in each non-trivial 
sector to quadratic order in the fluctuations [13, HI]. For a cylindrical shell it is essential to enforce the upper bound 
on the longitudinal momentum fraction < a << 1. Ignoring it and integrating without restriction over quadratic 
fluctuations about the stationary point at a = 0, the SCE of a cylindrical shell in fact vanishes [26l [3^ . However, 
restricting the quadratic fluctuations to < 1 essentially reproduces the held theoretic value for the Casimir energy 
of a metallic cylinder. The Casimir energy of a cylindrical metallic shell in a semi-classical sense is entirely due to 
Fresnel difltraction effects. The restriction of fluctuations to a < 1 in fact is required by causality: the spectral density 
otherwise is not analytic in the flrst quadrant. This violation of causality does not occur for any other fluctuation 
integral, which to leading order are evaluated without restriction. 

From a semiclassical point of view, the metallic cylindrical shell thus is a rather interesting geomery and conceptually 
more rewarding than the spherical one. Because a straightforward semi-classical evaluation to leading order resulted in 
a vanishing Casimir stress, it was previously believed_33] that this explained the exact vanishing of the Casimir stress 
on a dielectric-diamagnetic cylindrical shell (with equal speed of light on either side) to first order in the reflection 
coefficients [1, [13, [sl]. The non- vanishing Casimir energy of a metallic cylindrical shell was attributed to sub- leading 
diffraction effects. The last conjecture has now been verified, but the diffractive contribution due to classical orbits 
with just two reflections, the n = 1 contribution in Eq. (PT|) . is —0.0174076 .. .ficL/i?^. It does not vanish and 
is the largest contribution in magnitude, larger than the total self-stress of the metallic cylinder. Although we thus 
lack a semiclassical understanding of a weakly reflecting dielectric-diamagnetic cylindrical shell, it should be pointed 
out that the self-stress of a dielectric-diamagnetic spherical shell[6l| also is underestimated by a factor of 2.5 in this 
approximation. 

In the limit of very small reflection coeflicients it probably is impossible to ignore contributions from the exterior. 
It may be necessary to include diffractive contributions from paths that creep about the cylinder to semiclassically 
describe weakly reflective interfaces. Mathematically the exact cancelation in the cylindrical case is an addition 
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theorem for Bessel functions that expUcitly requires the exterior contributions [5a|. That perturbation in the reflection 
coefficients is quite deUcate becomes evident for dielectric cylinders [G^I and spheres [6^: to second order in the reflection 
coefhcients the results for the Casimir stress in this case are finite and comparable to the case where the speed of 
light is continuous across the boundary, but the total Casimir self-stress diverges logarithmically [H, if the speed 
of light in the interior and exterior do not match exactly. 
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